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ABSTRACT 

Detailed analysis of the coupled Dirac-Maxwell equations and 
the structure of their solutions is presented. Numerical solu- 
tions of the field equations in the case of spherical symmetry 
with negligible gravitational self-interaction reveal the exis- 
tence of families of solitons with electric field dominance that 
are completely determined by the observed charge and mass 
of the underlying particles. A soliton is found which has the 
charge and mass of the electron as well as a charge radius 
of 10~'^"'m. This is well within the present experimentally 
determined upper limit of ~ lO^^^m. Properties of these 
particles as well as possible extension to the work herein are 
discussed. 

1999 PACS numbers: 03.65. Ge, ll.lO.Lm, 12.20.Ds 

1 Introduction 

Through the years, a number of authors have attempted to avoid the problems inherent in the 
point-particle model by focussing upon finite soliton-like structures. Fields interacting non- 
linearly provide the binding without invoking any phenomenological elements. Einstein and 
Rosen pointed out many years ago that particles should be contained within a field theory 
and not exist as independent entities. Rosen Q made considerable progress in implementing 
such a program in a gauge-invariant manner by minimally coupling a scalar field to the 
Maxwell field. However, the soliton solutions yielded negative masses. Later 1^], neutral 
quantized particle states of positive mass were found and a more complicated model invoking 
up to three scalar fields coupled to the Maxwell field was shown to be capable of modeling 
the known massive leptons|Q]. However, the particles were spinless and the view then was 
that a subsequent quantization of the theory would induce spin. 

In 1991, one of the present authors[^ suggested an alternative route to elementary particle 
modelling, namely as solitons of Dirac-Maxwell theory. Since Dirac-Maxwell theory had 
been so successful in describing electron spin and magnetic moment, predicting the existence 
of the positron and refining the energy levels in interacting systems such as hydrogen, it 



seemed reasonable that this might successfully extend to a self-interacting soliton structure 
to model the elementary particles themselves. Spin would already exist in such a model via 
the spinor structure of the wave function. Shortly thereafter, such solitons were found and 
their properties studied[Q. A few years later, Lisi]^] independently discovered some of the 
results in|^]. Recently, there has been a revival of interest in this field and in particular, the 
issue of gravitational coupling in the Dirac-Maxwell system has been considered |^]. However, 
there was the misconception that gravitation was a necessary ingredient for the creation of 
the soliton. 

In this paper, we develop the essential results in[^ and@ and discuss the role of gravitation 
in soliton structure. The experimental inputs are the respective masses of the electron, muon 
and tau, their charge and as a constraint, the upper limit to their size which is ~ 10^^^ cm. 
The plan of the paper is as follows: in sec. ^, we set out the essential coupled Dirac-Maxwell 
equations to be solved. The structure of the Dirac wave function in spherical coordinates is 
given and particularized to the case of electric field dominance. The equation is separated 
in sec. ^ and we contrast the standard treatment in which a potential function is imposed 
such as in the case of hydrogen and the present case of the soliton where the derivation of 
the potential is part of the problem. The formal structure of the potential in terms of the 
Green's function is given. It is shown that there do exist spherically symmetric potentials for 
appropriate choices of quantum numbers. 

In sec. ^ the spherically symmetric energy-momentum tensor is derived. The relationship 
between the parameters in the Dirac equation and the physically measured quantities is 
discussed and the expression for the spatial spread of the soliton is given. The various 
constraints including singularity avoidance lead to the required boundary conditions for the 
problem. 

In sec. the results are presented. New variables of convenience for numerical integration 
are introduced. The parameters leading to twenty ground state solitons are listed. It is found 
that there is a critical range which leads to solitons within the experimentally observed upper 
limit to the size of the electron. Excited states are presented and the mass ratios are found. 

In the final section |^, the essential achievements as well as the limitations of the results 
are discussed. It is stressed that the solitons have been found without the requirement 
of significant gravitational interaction and it is conjectured that gravity will be significant 
for Dirac-Maxwell solitons when ejm ~ 1 in units for which G = c = 1. In cgs units, 
this is 2.58 x 10~^ esugm"^. By contrast, the e/m ratio for the electron is 2.04 x 10^^ or 
5.27 X 10^'' esugm^-*^ in cgs units. 



2 



2 Derivation of the Equations 

The field equations are obtained from tlie Lagrangian of quantum electrodynamics! 
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(1) 



where ip = ("^i, V'2i V'3> is the Dirac spinor, -ip = ip^j'^ = , -02 , — — V'Dj = ifi ^) 
is the electromagnetic four-vector potential and F^'^ = d^A'^ — is the Maxwell tensor. 

The 7^ are 4x4 Hermitian anticommuting matrices of the unit square 
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where / is the unit 2x2 matrix and the Ufc are the Pauli matrices 
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Variation with respect to A^ and if) respectively, yield the field equations 



0. 



(2) 
(3) 



If ip is chosen to be an energy eigenstate with energy E and one chooses a static charge 
distribution with a four-vector potential of the form 

A^' = [(l){r,e,v),A''{r,d,ip)), A; = 1,2, 3 

then the equations (0)-(|3|) are reduced to 

-ihca • V -I- a4^mc^ — ea-A + ecj) — Eip = (4) 

= -47reV'V (5) 
V X (V X A) = Airetp^aip (6) 

where = 'y^'j^. 

In spherical coordinates, {x,y,z) = (r sin 9 cos 99, r sin sin 93, r cos 0) , the Dirac wave func- 
tion has the structure llTC 
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where / = /(r), g = g[r) and the {Yl^{9, ¥?)}i,m is the set of orthonormal spherical harmonics 
defined for / = 0, 1, . . ., m = — /, — I + 1, ... A and 



YneM = y^fq:^ Pr(cos^)e--. (8) 

In addition, m is an integer such that — j < m + 1/2 < j; (m + 1/2)1% is the z-component of 
the total angular momentum. 

Consider the spinor with j = 1/2, / = and m = which implies from the above 
representation (0) 

47r^/;^Q;^/) = 2/(r)5(r) sin ^(— sin(/7, cos 0)"'". 
Resolving equations (^) and (^) into spherical coordinates gives 

VV = -e(/(r)2+5(r)2) 

V X (V X A) = 

f 

V X (V X A) , = 

V X (V X A) = 2e/(r)5(r)sin6'. 

Therefore, a four-vector potential of the form 

A^' = (0(r), -A^{r, 6) sin 99, A<^(r, 6I) cos 0) 
should be chosen where the components satisfy 

^ + 2^ = -e(/M^ + ,(,f) (9) 

ar^ r ar ^ ' 

Since the right hand side of equation (p!o|) is nonzero, the theory can only be exact if A^p is 
nonzero. However at this point we will impose the assumption of electric field dominance 
and hence the dominance of (j) over A or /(r) dominance over g{r). 

For the validity of the approximation A = 0, one radial component of the spinor must 
dominate over the other so that 

It will be demonstrated that such objects do exist within the non-linear field. With this ap- 
proximation the equations to solve reduce to a Dirac equation coupled to a Poisson equation: 

ificx • V + a^mc^ + V{r) \ ^ = Eil) (11) 
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(12) 



With these facts in mind, we now turn to the separation of the stationary Dirac equa- 
tion (PI) with respect to a general central potential and the derivation of the form of ^^^/^ 
for a general set of quantum numbers. 



3 Separation of the Equation 



The separation procedure follows that given in Bethe and Salpeter[|10|]. First one introduces 
quantum numbers I and j; I is the orbital angular momentum quantum number as well as 
being an integer > 0; j is the total angular momentum quantum number and can assume 
just the two values / + 1/2 and I — 1/2, (but only +1/2 for / = 0). The forms assumed by the 
four components of ijj are given explicitly in (^). 

The explicit form of the Dirac equation ( pT] ) for the four components of the wave function 

is: 
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(13) 
(14) 
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(16) 



Therefore, by inserting the assumed wave functions (0) into (|l3|)-(||) and using identities 
similar to ([A.6|) we find that the following two coupled equations between / and g hold: 
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where the new quantum number k is defined as 

f -I 



I 



1 for j = / + l/2 (/ = 0, 1,, 
for j = 1-1/2 (/ = 1,2,, 



(17) 
(18) 

(19) 

= V{r) and together 



These equations are valid for all spherically symmetric potentials V{r 
they replace expression (pT|). 

At this point, the standard procedure is to specify an external spherically symmetric 
potential, an example of which is the electrostatic potential energy of the proton-electron 
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interaction. That is, simply 



V{r) 



which is the fundamental solution of Laplace's equation||T 



(20) 



where 5^{r) is a three dimensional Dirac delta function centered at the origin. This is 
consistent with the far rang^ behaviour that we expect to find for the self-field of the fermion 
since, when we compare (pC|) with (12) we see that the fermion is treated as an object without 
structure through the equality, 

^P^^ = -6^{r). 

There is one additional problem that must be explored, namely how to couple relation 
^ to ([l7|)-(|l8|). This will be achieved in three parts. First, we find the Green's function for 
the equation (|5|). Second, we find an analytic form for the probability density V'^V' using (0). 
Once this equation is known, we can proceed to the third step which is to find V{r) by 
forming the convolution of the Green's function of step one, with the probability density of 
step two. 

The potential V{r) satisfies the Poisson equation (12) and by assuming that the solution 
is sufficiently regular, this can be converted to an integral equation p 



V{r) = -A-rre^ / G{r,r')'ip\r')'il;{r')dr' 



where G{r, r') is the Green's function of the Laplacian operator in three dimensions 



G(r,r') 
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(21) 



(22) 



With the Green's function determined, we can turn our attention to the probability den- 
sity. This is accomplished by using a pair of identities for the associated Legendre functions^ 

. 2 
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(23) 
(24) 



together with the definition of the spherical harmonics (|8|). The resulting expression for the 
charge density of the Dirac particle is given by 



21 + 1 



2 r 
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m] 



+ (l + m+l)\Y[' 



(25) 



^By far range, wo mean those distances much larger than the Bohr radius r ^ Ti^ /me^ . 
^Both Eqs. (p3|)-(p4|) foUow directly from Eqs. (8.5.1) and (8.5.3) of Abramowitz & Stegunp 
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when j = 1 + 1/2 and 
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(26) 



when j = I — 1/2. 

Therefore by using (^), (^) and (p5|), one obtains the expression 
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for the case j = I + 1/2. Similarly, with the use of (26), it can be shown that the potential 
V{r) takes the form 
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for the case j = 1 — 1/2. It is to be noted that the primed indices {l',m') correspond to the 
angular momentum of the particle, while the unprimed indices run over the complete set of 
permissible angular momentum quantum numbers. By performing the angular integration 
of the above formulae, one can immediately conclude that both of the above integrals vanish 
except when m = and / = 0, 2, . . . , 21' . This implies that 
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(27) 



where the cases j = 1 + 1/2 have been combined by the application of the definition of k' . 
Expression ( p7| ) replaces the equation (p!2|). When written in this form, it is clearly seen that 
the potential V{r) is not in general spherically symmetric. Table || lists the potential ( p7| ) 
for /' = 0, 1 and illustrates the fact that there exists spherically symmetric states with /' 7^ 0. 

A localized solution of this model must satisfy the field equations ( |T7| ) and ( [T^ ) for / and 
g and a given energy E where the potential is given by the expression (p7[). Moreover, it is 
required that the total probability 

(VIV) = E^V'.IV'.) = + r^dr < 00. 



i=l 
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Since the equations which describe the spatial evolution of the wave function ([I7|)-(p^) 
were derived under the assumption that the potential, V{r), is spherically symmetric, they 
are not valid for an extended Dirac particle in an arbitrary state of angular momentum. We 
have shown that there do exist certain choices of I and m where the probability density is 
spherically symmetric and it is these cases in which our primary interest lies. 

We can conclude that with the spinor representation given by (Q), there are essentially 



three differential equations to be solved simultaneously. Equations (|17D-(18) specify the 
spatial evolution of the wave function and equation (|2^) reflects the spatial extent of the 
self-field of the particle. A strategy for solving these intrinsically non-linear equations, as 
well as a few of their interesting properties, will be explored in the following sections. 



4 Boundary Conditions 



From the previous section we have found that the equations to be satisfied for a self-interacting 
fermion are equations (|T7|)-([T8|) and 
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where we have combined the j = I ±1/2 cases by using the definition of k. Since we have 



assumed that the potential V in equations (p^-(18) is spherically symmetric, this necessarily 
restricts the choice of / and m. Assume from this point on that / and m are chosen to satisfy 
this criterion. Consequently, equation (p8[) becomes 



{f + 9' 



(29) 



Since the soliton asires as a coupling between Dirac and Maxwell fields, the energy E that 
appears in the Dirac equation is not the total energy of the particle. The total energy can be 
obtained by calculating the Tq component of the energy-momentum tensor. For our field, the 
Lagrangian is given by equation (^ where A^^ is the vector potential of the electromagnetic 
field. One generates the symmetric energy-momentum tensor directly from the Lagrangian 
in the form [ 14 1 

(30) 
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Applying (|0D to ^ yields 
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Further simplification gives 
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This yields an expression for the total energy, -Etot; of 
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(31) 



where dVoi is an infinitesimal volume element. This total energy should be associated with 
the observed mass of the particle as Etot = rnc? . There is still sufficient freedom remaining 
to set limr^oo V{r) = because the spinor is invariant under the transformation 

V^V + 13; E^E + /3 

for any real- valued (3. 

The mass m and the charge e that appear in the Dirac equation are not necessarily the 
experimentally measured quantities just as the charge that appears at a vertex of a Feynman 
graph is not the experimentally measured charge of the particle. Because of this, we will 
replace the m in (p!7|)-(p^) by the symbol /U. In addition, the e in ( p9| ) will be replaced by 
an e. The symbols m and e will be reserved for the physically observed quantities. With 
these substitutions, we convert to a set of variables whereby equations (p!7[), (|l8|) and ( pl| ) 
are independent of any physical constants. The particular transformation chosen is 



f = VE, g = r]G, r = —, 

lie 



E = Xnc^, V 



where r/^ = ic'c^ / e^fi? . These redefined variables have the following dimensions in terms of 
length (L): 
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This yields the transformed equations: 
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(32) 

(33) 
(34) 



where V'^ is now the Laplacian with respect to the x coordinate. The mass of the soliton 
comes from the transformed version of the total energy expression ( |3ll) , 
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and the total charge is given as the integral of the charge density 

he roo 



= ejpdVoi = ^j^ (f^ + G') dx. 



(36) 



We will show that if the charge e is replaced by e, that the / component of the spinor 
greatly dominates the g component. By substituting e = e and choosing a value for m, the 
value of fi can be determined numerically once the spatial extent of the soliton is known. In 
this case, the expectation value of the radius of the particle becomes 



^^(^^^^ V^Ux^dx 



V^Ux^dx 



I 



V'^Ux^ dx 



To stay within the current experimental bounds of the mean charge radius, this value must 
be less than Vexp which is < 10~^^m in the case of an electron. Hence, 



/ vWcia;<^-^ / VWdx 

Jo Te rUe Uo 



where Vg is the classical electron radius re = e^/mgC^. 

Since we know that U has zero slope at x = and that it must behave as N/x for large 
argument {N is the amount of enclosed charge), we assume as a first approximation, that 
U can be represented as the electrostatic potential produced by a sphere of radius Rq with 
uniform charge density. Therefore, 



U{r) 



N 
N 



111 
2Ei 



for r < Ro 
for r > Rq. 



(37) 



With this representation, one finds that 

dretneRo 
(r) = 

which means that since < (r) < rexp, we can conclude that 



< 



^0 



< 



4 r, 



exp 



3.088 X 10-^ cVMeV 



in the case of the electron. 

Let Rq be defined as the effective range of the non-Coulombic behaviour of the potential 
energy so that for x > Rq, U{x) N/x. Since U is & solution to a Poisson equation with 
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a negative definite charge density, U{0) must be larger than U{Ro). This can be quickly 
verified by considering the opposite. If U{0) < U{Ro) then there exists some r G (0, i?o) such 
that U'{r) = 0. Therefore integrating from to r, one obtains 



r^U'ir) = = - (^F^ + G 



dx 



which is clearly a contradiction. 

To determine the initial values of F and G, one simply eliminates either F or G from (I 



33), say F, which leads to a second order equation for the other, namely, 

G" + PG' + QG = 

where both P and Q are functions of U, U', k, A and x. To avoid a singularity in the 
potential U{x), it must be both bounded and have zero slope in a neighbourhood of the 
origin. Moreover, both F and G are bounded in this same neighbourhood. From this, it is 
easy to verify that 

' K = -l 

F(0) = { arbitrary k = +1 (38) 
V other k, 

arbitrary k = —1 
G{0) = { K = +l (39) 

V other n. 

Furthermore, by examining the indicial equation, it can be shown that no fractional powers 
exist in a power series solution of either F or G about the origin x = 0. 
Summarizing the boundary conditions: 

U{x) < U{0) < oo, xe [0,oo), 

U'{0) = 0, 

together with the conditions (|38|), (|39|). For the case k = —1, the initial values of U, G and 
the energy A are determined by the requirement that the wave function tp, and hence both 
F and G, vanish exponentially as a; ^ oo. 



5 Results 

In the search for numerical solutions it was specified that k = —1 and A = 1 giving the set of 
differential equations 

^ = [2-U{x)]F{x) 
dx 

dF 2 

^ = F{x) + Uix)Gix) 

dx X 
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To find a soliton, the values of F{0), G(0) are specified and a searcli is made for the value 
of U{0) whereby lima;_^oo xF{x) = and lima;_+oo xG{x) = 0. Only values of G(0) > 
are considered because the equations are symmetric under the transformation G — > — G, 
F — > —F, U ^ U. In a neighbourhood of a ground state soliton, the radial probability 
density is numerically seen to have a single well-defined minimum value for x > 0. The 
choice of K = —1 gives the initial condition F{0) = 0. 

The choice of A = 1 is simply a numerical convenience. Outside the neighbourhood of 
a soliton it is expected that the potential will behave as U{x) ~ A + B/x for large x. The 
value of A should have been chosen so that the asymptotic behaviour of the potential U{x) 
is purely Coulombic in nature. By defining a shifted potential Uix') = Uix') — lima;_j.oo 
this value of A must satisfy 1 — U{x) = A — U{x). Therefore after a soliton is found the value 
of A is given as A = 1 — lim^^oo ^^(2;). In addition, the starting value of U{x) is given by 
U{0) = c/(o) + A - 1. 

Using the redefined value of A the observed charge and mass of the particle are compared 
to the values used in the Lagrangian by using the expressions (^) and (^) respectively. By 
defining 

r = (f"^ + G'^)x^ dx, X = (f"^ + G"^) x^ dx, 
the charge ratio e/e is given as 

e a 

where a is the fine structure constant. The mass ratio fj,/m = V'^/aS and the expectation 
value for the radius of the soliton is 

_ /(/^ + a'^y dr _ n J{F^ + G^)x^ dx _ fmA £X 

Both of the quantities V and X are positive. However, depending upon the value of A, £ 
could be positive, negative or even zero if the electromagnetic and "bare mass" terms in the 
energy exactly cancel. A negative value for £ will give an unphysical negative value for the 
observed radius (r). Because of this ambiguity, both the value of (r) and the particle width 
Ar = V {r''^) — {r)'^ are presented. Tables ^ and |^ respectively list the numerical parameters 
and the observed properties of a number of ground state particles found where m was taken 
to be the observed mass of the electron mg. 

Figure |l] illustrates the radial behaviour of F and G for the case e = e (i = 1). It is to 
be noted that for x > 0, -F is much larger than G and as a consequence, FG ^ + G^. 
In fact, G is so small that it resembles a straight line along the x axis. This supports the 
argument that the four-vector potential can be reasonably approximated with only a radial 
A? component. 
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The characteristics of a typical sohton with e 7^ e is illustrated with the choice e/e = 
454.8 (i = 19). In this case the potential plays a much more dominant role in holding the 
particle together than in the case e = e. However, since in this case the approximation of 
FG <^ F^ + is violated, one would have to solve the full model (equations (pD-(p^)) to 
properly analyse this situation. This would be a far more complicated problem. Figure |2| 
illustrates the radial components of this spinor and it shows that the magnitude of G is 
now comparable to the magnitude of F. Table ^ also shows that the choice of e = 389. Oe, 
fi = 2.360 X lO^^me {i = 15) yields a soliton with an expectation value for the radius of 
5.05 X 10^23^^ 'pj^jg 

size is well within the present experimentally determined upper limit 
for the electron radius of ~ lO^^^m. 

These equations also exhibit excited states. The n*^ excited state of our field is charac- 
terized through the functions Fn{x), Gn{x) and Un{x) for which the Gn component crosses 
the abscissa n + 1 times while the Fn component crosses it n times. Once the ground state 
solution is found, the value of /i can be determined through equation (36). The corresponding 
v}^ excited state is that excited state with the same observed charge ratio, e/e, as the ground 
state. Therefore, in this interpretation of the theory, the ratio of the mass of the n*^ excited 
state to the ground state is given by the expression 

Figure 1^ shows radial probability density of the first three states for the case G(0) = 1. Each 
of these solitons has a different value of e/e. 

Figure § illustrates the behaviour of the mass ratio, /i/m, as a function of the charge ratio 
e/e for the ground state and the first two excited states. For each class of particles there is 
a charge ratio where the electromagnetic and bare mass components of the energy balance 
making f = 0. At this value of e/e, the mass ratio /i/m — > cxd. At charge ratios less than 
this critical value the mass ratio is negative whereas charge ratios above this critical value 
result in a positive value of ix/m. There is numerical evidence that each class of particles 
has an upper bound for the charge ratio. Above this maximum charge ratio we were unable 
to find any solutions such that lim^^^oo xF[x) ^ or lima;^oo xG{x) — > 0. This necessarily 
restricts the definition of the mass ratio defined above. Figure ^ also illustrates the fact that 
at moderate charge ratios, the electromagnetic field does not contain an appreciable amount 
of the particle energy resulting in the behaviour \^/m\ ~ e/e. 

The mass ratios of the first and second excited states with respect to the ground state 
solutions are shown in figure ^. This ratio is only defined up to a maximum value of e/e since 
beyond e/e ~ 550, a ground state fails to exist. For excited states, this maximum admissible 
charge ratio increases. This implies that for a fixed value of e/e there may not exist a ground 
state solution, but there will be arbitrarily many excited states. As is readily apparent from 
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figure ^, the only appreciable mass splitting occurs for large charge ratios. However, it is 
precisely for large charge ratios where our approximation that FG ^ + breaks down. 

6 Concluding Remarks 

We have seen that spherically symmetric Dirac-Maxwell solitons can be constructed and with 
a charge and mass to model the electron successfully. However, it should be noted that the 
higher energy excited states of this form did not yield the large mass separations of the muon 
and tau relative to the electron in this model. The search thus far has been restricted to 
spherical solitons. It is conceivable that a relaxation of this restriction or some other change 
in conditions would increase the mass splitting. In any event, we have shown that Dirac- 
Maxwell solitons exist and are capable of modelling an electron where the charge-to-mass 
ratio is the observed ~ 10^^ in units in which G = c = 1. Furthermore, we have found a 
charge-to-mass ratio that simultaneously yields the observed charge and mass of the electron 
as well as exhibiting a degreee of compactification that is well within the current experimental 
upper limit. Finster et al.[^ have considered Einstein-Dirac-Maxwell (EDM) solitons and 
concluded that it is the interaction with gravitation which is responsible for the existence 
of bound states. However, we see here that bound states exist with negligible gravitational 
interaction. While the e/m ratio at which significant gravitational coupling sets in is yet to 
be determined for EDM solitons, it is our conjecture that this will be so at the same level that 
was found earlier in the case of minimally coupled scalar interaction namely for e/m ~ 1. 
The known fundamental charged particles of nature, on the other hand have enormous e/m 
ratios. 
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A Derivatives of f {r)Yi^ {6 , (p) 

In the Dirac wave equation, all of the derivatives are with respect to Cartesian coordinates. 
We can change to a spherical polar representation via the transformation 

X = r sin 6 cos (/? 
y = r sin 6 sin ip 
z = r cos 0. 

By applying the chain rule, it is trivial to show that this changes the first order partial 
derivatives via 

^^"'^ ^ (A.1) 
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= sin 6 


cos (f 


d 
dr 


cos 


9 cos if d 


dx 




r do 
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= sin 6 


sin (p 


d 
dr 


cos ( 


9 sin 99 d 


dy 




r do 
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dz 


= cos t 


dr 


sin 6 
r 
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89' 





r sin 9 d(p 



cos ip d , , s 

r sm t> dip 

(A.3) 

If the functions ipj{j = 1,...,4) from expression (^ are substituted into (p^-(16), and if 
one uses the formulas given in Bethe and Salpeter [l^ for the derivatives of a function of the 
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form f{r)Y[^{9, (f) with respect to x, y, and z, one finds a coupled pair of first order ordinary 
equations for /(r) and g{r). 

For example, in order to calculate 



we first require the identities 



cos ^P;™ (cos ( 



21 + 



- [(/ - m + l)Pi'^^i(cos e) + {l + m)Pl^^ (cos 9)] 



(cos 9) = ^ [/(/ -m + l)P,'^i(cose) - (/ + 1)(/ + m)P,-i(cos0)] 



which can both be verified through the use of Rodrigues' formula 



( 1 \m ^il+m 

,.'2\m/2 " {,.2 



Writing Y[^{9, ip) as a function of P[^^ by using (^ gives the relationship 



(A.4) 
(A.5) 



^Jf{r)Yr{9,^)] 



I2l + l{l-m)\ 
Att {l + m)\ 



cos9 P[^{cos9)^ - sm9 4;Pr{cos9)- 
dr d9 r 



By substituting (A.4- A.5) in the above, collecting terms, and applying the definition of 
Yi^{9,Lp) once again, one obtains the simplification 



^[f{r)Yr{9,^)] 
oz 



^{l-m + l){l + m + l] 
(2/ + l)(2/ + 3) 



+1 



' (/ - m){l + m) 
{21 - l){2l + 1) 



dr r 



(A.6) 



Similar relationships for ^ it i-^ can be found in Bethe and Salpeter^, but there is a very 
elegant way to derive these operators by applying the Wigner-Eckart theorem. 

First, we evaluate the matrix element {I 0\Vo\l 0) of the gradient operator, which is an 
example of a vector operator. Specifically, 



Vn 



d_ 

dz 



, V± = T 



1 / d 



±i- 



d 



a/2 \dx dy J 



Since 

Vo/(r)y,° 



/ + 1 



V(2/ + 1)(2/ + 3) 



dl_l_ 
dr r 



f 



+ 



I 



y/{2l-l){2l + l) 



-yP 



dr r 



^See formula (A. 38) and (A. 39) respectively in Bethe and Salpeter. 
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for the special case of (A. 6) where m = 0, we have 



{I' 0\Vo\l 0) 



v/(2/ + l)(2/ + 3) 



dr r 



I 



v/(2/-l)(2/ + l) 



dr r 



>l-i- 



Now, we are at a point where we can use the Wigner-Eckart theorem. By inspection, the 
general matrix element is given by 



{l'm'\V^\lm) = i-iy-' 



i-iy 



I' 1 I 

-m' fi m 



{l'\\V\\l) 



I' 1 / 

-m' fj, m 



I' 1 / 




-(/' 0|Vo|/ 0). 



After evaluating the 3 — j symbols, one can quickly verify the following equations. 
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(?-m + l)(/-m + 2) ^^^_i^^ 
(2/ + l)(2/ + 3) '+1 ^ '^^ 

4f / + 1 



(A.7) 



(A.8) 



---/ 
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Linear combinations of (A.8) and (|A.9D yield the derivatives with respect to x and y. 
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Stales \l',m',j') Con'(;spondiiig Polculial V{r) 


|0, 0, 1/2), |0, -1,1/2) 
|1, 0, 1/2), |1, -1,1/2) 

|l,l,3/2),|0,-2,3/2) e2 
|1, 0, 3/2), |0, -1,3/2) e2 


/o-i(3cos2 0'-l)/2 

/0 + |(3COS2 e' -1)I2 





Tabic 1: The Dirac-Maxwell particle self-field potential. 
The self-field potential energy for a Dirac-Maxwell particle in the states ^ = 0, 1 where 
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Table 2: Numerical parameters for a set of various ground state particles. 
For each particle, the value of G{0) is selected and one searches for the value of U{Q) + X that gives 
a bounded solution. The physical parameters are computed from the solution defined on x G [0,a;max]- 
G*(0) =0.4273589430. 
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Table 3: Corresponding observable quantities for a set of various ground state particles. 
The values (r) and Ar are measured in meters and are computed from a soliton defined onx G [0, Xj^ax]- 
For this calculation it is assumed that m = mg. 
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X 

Figure 1: F and G components of the wave function for the case e/e = 1. 

Shown here is the radial dependence of the F and G components of the soliton. Note that G is much 
smaller than F. G is barely discernible above the x axis. 
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Figure 2: F and G components of the wave function for the case e/e = 454.8. 

Shown here is the radial dependence of the F and G components of the soliton. The magnitudes of F 
and G are now comparable in contrast to the case when e = e. 
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Figure 3: Excited states of the theory. 

This figure shows the radial probability density of the first three particle state in the case G(0) = 1. 
These particles have different charge ratios e/e. There exist excited states beyond the ones illustrated. 
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Figure 4: Dependence of the mass ratio as a function of the charge ratio. 

Shown is the dependence of the mass ratio fi/m as a function of the charge ratio e/e for the ground 
state and first two excited state solitons. For each class of particles, there is a maximum charge ratio 
beyond which no solitons were found. 
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Figure 5: Mass ratios of the first and second excited states with respect to the ground state. 

There are essentially two regions of interest. For moderate charge ratios, the value ofmn/ruo ~ 1 
with the mass ratio of the excited state n = 2 slightly larger than for the n = 1 state. Beyond the point 
where the ground state mass ratio becomes unbounded, the mass ratios begin to split. In this region, 
the \mi/mo\ ratio exceeds the |to2/too| ratio. In this region the approximation FG F^ + is no 
longer valid. 
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